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SUMMARY 


A theoretical investigation is undertaken to study the buckling of a 
thin, rectangular, homogeneous, isotropic flat plate augmented with 

a centrally located longitudinal stiffener. The influence on overall 
plate efficiency of the use of advanced composite materials in the 
stiffener is developed and discussed. Results of the analysis are 
presented in the form of charts to demonstrate the effects of fibrous- 


composite reinforcement of the stiffener on the critical-stress-to- 


weight ratio of the plate. 
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ListeOF SYMBOLS 
arbitrary constant in lateral deflection func- 
Elon. ii. 


cross-sectional area of longitudinal stiffener, 
in.© 


width of flat plate, in. 


flexural stiffness of homogeneous, isotropic 
plate, lb-in. 


flexural stiffnesses of orthotropic plate in 
axial and chordwise directions, respectively, 
Tbereire 


torsional stiffness of orthotropic plate, lb-in. 


Young's modulus of homogeneous isotropic plate 
and stiffener, Dy ar 


extensional stiffenesses of orthotropic plate in 
axial and chordwise directions, respectively, 


iby an. 


effective flexural stiffness of composite-reinforced 
stiffener, lb-in.¢ 


shear modulus of homogeneous isotropic plate and 
stiffener, lb/in.2 


shear stiffness of orthotropic plate, lb/in. 


effective torsional stiffness of composite-reinforced 
stiffener, 1b/in.© 


stiffener height, in. 
moment of inertia, in. 
torsional stiffness constant, in.4 


arbitrary nondimensional integer 


buckling coefficients in axial and chordwise direc- 
Elomcr becpeceivery 


length of flat plate and stiffener, in. 
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plate bending moments in x- and y-directions, 
respectively, lb-in. 


number of buckling halfwaves in x- and y-direc- 
tions, respectively 


middle surface stress resultants, 1b/in. 
applied load resultants, 1b/in. 
arbitrary nondimensional integer 

torque, lb-in. 

plate thickness, in. 


displacements of point on middle surface of plate 
in x-, y-, and Z-directions, respectively, in. 


total strain energy of plate and stiffener, lb-in. 
strain energy of plate, lb-in. 

strain energy of stiffener, lb-in. 

potential of applied loads, lb-in. 

Wadith OL “actwal) stiiienen, in. 


equivalent width of composite portion of stiffener, 
in. 


plate coordinates, in. 
buckle aspect ratio of unstiffened plate 
beam twisting strain, in./in. 


transverse shearing strains in x- and y-directions, 
respectively, in./in. 


plate shearing strain, in./in. 
burac idelllta Gumetvonsay ane 
Kronecker delta function 


membrane strains in x- and y-directions, respectively, 
Was bide 


angle of twist, radians 
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SUBSCRIPTS 


base 


comp 


Poisson's ratio for homogeneous isotropic plate 


Poisson's ratios for orthotropic plate in x- and 
y-directions, respectively 


density, 1b/in.? 
average density of plate and stiffener, iy en 
a 


stiffener middle-surface stress, 1b/in. 


plate membrane stresses in x- and y-directions, 
respectively, 1b/in.© 

er eS 

critical stress of plate, 1b/in. 


stiffener shear stress, 1b/ in.” 


pertaining to base material (aluminum) 


pertaining to composite material (boron epoxy or 
graphite epoxy) 
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INTRODUCTION 


The trend to higher performance aircraft in all speed regimes has led 
designers to seek methods of decreasing the structural weight of air- 
craft while maintaining or increasing the stress levels on the mate- 
rials involved. One of the most promising means of accomplishing this 
objective has been the introduction of advanced composite materials to 
either augment or replace conventional structural materials (for exam- 
ple, aluminum). One type of composite material consists of the in- 
corporation of high-strength, high-modulus filaments into a compatible 
matrix. Several material system classes currently of active interest 
include graphite epoxy, boron epoxy, and metal matrices having metal 
fibers. As pointed out in Reference 1, design and structures engi- 
neers, using composites, have the ability to build directional prop- 
erties into the structure which, in turn, results in improved efficiency. 
Weight saved, with its direct effect on the design index of strength/ 
weight, is the most attractive feature of the use of high-performance 
composites. It has not been uncommon for structures-materials engi- 
neers to predict weight reductions of as much as forty percent in total 
systems weight through the potential use of advanced composites in both 


primary and secondary structures. 


With advanced filamentary composites, increased design flexibility of 
significant degree exists that is not possible with conventional metal 
structures. However, as with many new developments, work in the com- 
posite field has been primarily empirical, while theory and analysis 
have lagged somewhat behind. Despite the potential of composites, 
evidence of their use appears only in the reasonably small structural 
elements that can be designed by test (for example, stabilizers, wing 
tips, slats, and rotors). The trial and error procedure, especially 
with materials of such high performance possibilities as composites, 
can lead to excessive costs due to both inherent high material costs 
and eventual appearance of structural behavior neither anticipated nor 


properly estimated. One such case is presented by Durlofsky and Mayers 





in Reference 2. Another such case exists with the use of advanced 
composites to enhance the buckling efficiency of stiffened plates 
utilized in the construction of wing, fuselage, and empennage struc- 
tures of aircraft. Haereeene studied the effect of an all-composite 
stiffener in the reinforcement of a simply supported aluminum plate 
without considering the possibility of an optimized composite-metal- 
stiffener combination. His results indicate that some caution must 
be exercised in using high-performance composites in combination with 


conventional structural concepts and prediction methods. 


Earlier studies of conventional stiffened plates under compression, 


p) 


such as those of parece Seide and Stein,” and Gomza and Stein, 
neglect the effects of torsional stiffness of the stiffener in es- 
tablishing buckling criteria. The justification offered is usually 
that the torsional stiffness of the stiffeners is relatively small 
when compared with the bending stiffness and that vanishingly small 
torsional stiffness of the stiffeners is consistent with the appear- 
ance of buckle nodes at stiffener locations in the chordwise direc- 
Cton. Such justification, of course, is appropriate for the discreve 
stiffener spacing of low subsonic speed flight vehicles. For the 
higher speed ranges, stiffeners are more closely spaced and the 
stiffener torsional stiffness becomes significant in buckling anal- 
ysis. “Thus, the introduction of composite materials into all or pare 
of a stiffener can possess ramifications not evident in the use of 
conventional materials. For example, although Young's modulus OD) 
for a composite may be as much as four times greater than that of a 
structural metal, the shear modulus (G) may be lower by approxi- 
mately the same factor. Therefore, neglecting torsional rigidity 

in a composite-metal system with closely spaced stiffeners may in- 


troduce unsatisfactory errors in buckling calculations. 


It is the purpose of this investigation to (1) present the theory and 
an approximate analysis for determining the stability of a simply sup- 


ported, flat, rectangular orthotropic plate with a centrally located, 
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composite-augmented, longitudinal stiffener and (2) determine for a 
square plate the effects of varying the percentage (by volume) of 
composite material in the stiffener with stiffener torsional rigidity 


included. 





GENERAL THEORY 


STRAIN-DIS PLACEMENT RELATIONS 


The strain-displacement relations for small strain, large displacement, 
and moderate rotation behavior of an orthotropic plate (see Figure 1) 


can be obtained from Reference 7 as 


it 
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With a single stiffener buckling into no more than 2 half-wavelengths 
and a homogeneous, isotropic base plate, however, interlaminar and 
transverse shear effects become of higher-order significance and they 
are assumed to be negligible throughout the analysis. Thus, the trans- 
verse shear angles se and Ty are set equal to zero. Further, since 
buckling is of interest, the stretching of the middle surface of the 


plate (and stiffener) may be neglected; thus 
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x » XX 
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STRESS-STRAIN RELATIONS 


For an orthotropic plate in which the stresses do not exceed the elas- 


tic limit, Hooke's law is applicable; it is expressed in the following 
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form relating normal force per unit length of cross section to strain: 
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A single longitudinal stiffener possesses the same stress-strain re- 


lationships as a one-dimensional beam; that is 


G “= “Ee 
x 


(4) 


POTENTIAL ENERGY 


The strain energy stored in an elastic plate of constant thickness due 


to bending and twisting may be written as 


aS & eee 
U = — Nc + Gat N dxdydz 
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Substitution of Equations (2) and (4) into Equation (5) and subsequent 


integration over the thickness of the plate yields 
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D 
+ — 2 fF +on we | dxay (6) 
1-up TY XY ,xXY 
x y 
With the use of the Dirac delta function, the potential energy of a 
single centrally located symmetric longitudinal stiffener (see Fig- 


ure 1) may be described over the same interval as 
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Substitution of Equations (2) and (4) in Equation (7) and subsequent 


integration over the stiffener height yields 
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If an axial, inplane load, positive in compression, is applied to the 
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dxdy (8) 
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edges of the plate, the potential of this applied load (including the 


end load on the stiffener) is 
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However, since in the absence of middle surface stretching (<. =€ = 


Ves = 0 when Z = 0), 1t can be observed trom Equation (1) that 
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Substitution of Equation (10) in Equation (9) yields 
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Finally, the total potential is the sum of Equations (6), (8), and (11) 


and becomes 
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VARIATIONAL PRINCIPLE AND GOVERNING EQUATIONS 


In accordance with the potential energy principle, the total potential, 
U+V_, must have a stationary value with respect to admissible var- 
iations in the lateral deflection, w , in order for equilibrium to 


be achieved. This variation leads to the equilibrium or Euler equation 
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of the variation, plus the applicable boundary conditions. The var- 


iation of Equation (12) leads to the following equilibrium equation: 
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The associated boundary conditions for the case of single support are 


obtained as 


D 
x 
w= 0 and - W =M = 0 @ x = O,L 
l-up XX My 5) x ‘ 
x y 
(14) 
D 
w= 0 and - —~— lw + Uw =M = 0 @ y = 0,b 
eee Oy, x) ook y 
x y 


The solution of Equation (13), consistent with the boundary conditions, 
Equation (14) leads to the deflection function, w . This, however, 
is a difficult procedure in general and is not the approach followed 

in this investigation. Rather, the direct method of the calculus of 
variations (Rayleigh-Ritz method) is used to determine the buckling 


criteria for w expressed in terms of a series of arbitrary parameters. 





METHOD OF SOLUTION 


With the total potential energy given by Equation (12) in terms of the 
plate deflection w , the direct method of the calculus of variations 
(Rayleigh-Ritz method) is used to determine the minimum buckling load 


for varying stiffener size and composite content. 
GENERAL CASE 


A deflection function satisfying the boundary conditions for a flat, 


simply supported, stiffened rectangular plate under axial compression 


is 
- > _. TX .. gory . 
w = Ae sina sink. (15) 
n 
Substitution of Equation (15) into Equation (12) (with Ny = Se =O) 


and subsequent integration over the surface of the plate yields (see 


Appendix I) 
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The vanishing of the first variation of the total potential energy with 
respect to An , the nondimensional arbitrary displacement parameter, 
yields the equilibrium condition and a criterion for obtaining the 


buckling load. This procedure, as shown in Appendix I, yields 
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HOMOGENEOUS , ISOTROPIC SQUARE PLATE 


A homogeneous, isotropic square plate with an equivalent stiffener (to 
account for varying amounts of composite material) has been selected 
as the specific case for examination in the present investigation. A 
square plate with one centrally located longitudinal stiffener under 
these loading conditions would buckle in one of two configurations: 
(1) one half-wave in each direction (m=ne=1) or, (2) two half- 
waves in each direction (m=n= 2) . Any other combinations of 

m and n would not lead to the minimum value of the buckling co- 
efficient k, . As noted from Equation (17) the bending stiffness 
term (EL)/Db of the stiffener is effective only during the 1,1 
mode, and conversely, the torsional stiffness term (GJ) /Db 1S 7ef- 


Bective only during the 2,2 mode. 


For an approximate solution, the first term of the deflection function, 
Equation (15), could be used. However, this would allow only for a pure 
sinusoidal half-wave deflection pattern of the plate in the chordwise 
direction. As stiffener size or stiffness is increased, rotation at the 
inflection point of the 2,2 mode is somewhat restrained, indicating the 
need for at least a two-term deflection function; that is, n= 1,3 for 
the 1,1 mode and n= 2,4 for the 2,2 mode of Equation (17). The 


result of carrying out this procedure (see Appendix I) is 
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for the 2,2 mode, where D is the bending rigidity for a homo- 
geneous, isotropic plate. The quantities (EI) and (GJ) are the 
bending and torsional stiffnesses, respectively, of the equivalent 
longitudinal stiffener. The procedure for determining these values 
for varying amounts of composite material in the overall stiffener 


is given in Appendix II. 


The minimum values of k.. may be calculated from Equations (18) and 


(19) in each of the two modes of the square, stiffened plate for 





varying values of (EI)/Db , (GJ)/Db , and (A_)/bt ihe ier ered 
stress, O. ; is then calculated from the relationship 
k nD 
o = (20) 
Cie b“t 


The average density of the plate-stiffener system is a nondimensional 
measure of the relative weight of the combination with the introduc- 
tion of composite materials into the system. This factor is calculated 


by averaging the densities of the system as follows 


Pave = Pcomp (# ComP) + Prage (@ base) (ey 


where (% comp) is the percent by volume of composite in the entire 
plate stiffener system, and (% base) is the percent by volume of the 


Base metal used (aluminum, steel, etc. ). 


The buckling efficiency of the stiffened plate combination (Teed 

is plotted versus (A_)/bt for varying amounts (by volume) of composite 
material in the stiffener in Figures 2 through }. Aluminum is the base 
material, with either boron epoxy or graphite epoxy used as the advanced 


composite material. 


Beene determined that the optimum longitudinal stiffener for a flat 


plate has the relative proportions (1) stiffener width equal to twice 


Te 





plate thickness, and (2) stiffener height equal to eight times the 
plate thickness. These ratios are utilized as starting values for 


the present investigation. 


1) 








RESULTS AND DISCUSSION 


An investigation of the axial compression buckling of a simply 
supported flat plate with a central longitudinal stiffener has been 
undertaken to estimate the effects of composite reinforcement of the 
stiffener with stiffener torsional stiffness effects included. The 
minimum or critical values of k_ , using Equations (18) and (19), are 
calculated for the two modes of a stiffened, square, homogeneous, iso- 
tropic plate for varying values of (EL) /Db and plotted in Figure 5. 
These values are in excellent agreement with those of Seide and Stein? 
ror ‘the’ m= 1 , n = case, but the (Ke) values of this investi- 
gation are slightly higher than those of Reference 5 for the m= 2 , 
n= 2 case. This discrepancy is due to the inclusion of the torsional 
stiffness of the stiffener in the present analysis. Without torsional 
stiffness effects, (Kx) for the m= 2 ,n=2 case remains con- 
stant at a value of 16 despite the size of the stiffener. Timoshenko 
and cous” in the case of a plate simply supported on three edges with 
one of the unloaded edges clamped under compressive loading show the 
maximum value of (Ky mE to be 22.6. This would be analogous to the 
present problem when the torsional stiffness is infinitely large. The 
two-term deflection function used in the present investigation will 
result in values of (k _ higher than that arrived at by Reference 8, 
but only when the stiffener size is excessively large and out of the 
area of practical interest. Equations (18) and (19) give values of 
(kx) which are sufficiently accurate, as the primary aim of the 
present report is to investigate the trends of plate efficiency rather 


than to calculate actual design data concerning buckling loads. 


The primary quantity for demonstrating plate buckling efficiency in this 
investigation is Oia Dinan » where Tor is the critical stress of the 
plate determined from Equation (20), and Qe is the average density 
determined from Equation (21). Figure 2 is a graph of oo 1/ Pave for a 
square aluminum plate with a centrally located longitudinal stiffener 


augmented by varying amounts of graphite-epoxy composite. The buckling 
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efficiency has been calculated for varying stiffener size. The sharp 
knee in each of the curves results from the transition of the plate 
buch ng’ pattern from a “m= 1, n= 1 configuration to a m=e2 , 

n= 2 mode. To the left of this knee, or in the 1,1 mode, sig- 
nificant increases in buckling efficiency result from the replacement 
of aluminum with graphite epoxy in the stiffener. This is due to the 
increased flexural stiffness of the stiffener, since the graphite epoxy 
has a greater value of Young's modulus than does aluminum. However, to 
the right side of the knee, in the 2,2 buckling configuration, the 
plate buckling efficiency is actually reduced with the addition of 

the graphite epoxy. Once the 2,2 condition is reached, the stif- 
fener lies at a node of the plate, and the flexural stiffness of the 
stiffener has no effect on the magnitude of the critical stress. Only 
the torsional stiffness of the stiffener is a factor in this range and, 
Since the shear modulus of aluminum is approximately three times that 
of graphite epoxy, the overall buckling efficiency decreases with the 
addition of the composite material. As an example, a plate with an 
all-aluminum stiffener and a value of (A_)/bt of 0.205 attains a 
buckling efficiency parameter value of 41.0 x 10" im: (Beant). |e 
the aluminum stiffener is then replaced by an all composite stiffener, 


this parameter value is decreased to 38.8 x 10° in. (Point B). 


It may also be noted from Figure @ that in the 1,1 configuration, as 
the percentage of composite in the stiffener increases, the percent 
increase in the buckling efficiency decreases. For a plate with a 
ie. 

value of (A_)/bt of 0.08 , o,/p,,, equals 15.7x 10 in. for 
an all aluminum stiffener. The replacement of 10 percent of the stif- 
fener cross-section with graphite epoxy raises the value of oa /e 

i cr’ "ave 
to 19.5 x 10 in., a percentage increase of almost 25 percent. Raising 
the percentage of composite to 50 percent gives a buckling efficiency 
parameter value of 28.3 x 10 in. , or an increase of 80 percent 
over the all aluminum stiffener. An all composite stiffener has a 
value of 431.6 x 10" , Or a@ littlesgreater thanga 1COgpen- 


o 
es 
cent increase over the metal stiffener. In another case, where the 
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size of the stiffener is increased to give a value of (A_)/bt of 
O.14 , a 10-percent-composite stiffener has approximately the same 
value of Ty! Pave as Che all-composite stiffener, although the 


latter structure will be lighter in weight and much higher in cost. 


As mentioned previously, past studies having to do with the buckling 
properties of all-metal stiffened plates have generally neglected the 
torsional rigidity of the stiffener. Figure 3 is a graph of the same 
plate-stiffener combination as that of Figure 2; however, the shear 
modulus of graphite epoxy is assumed to be identical to that of alu- 
minum (G = 4.0~x LO veya. =) . In this case, the replacement of a 
percentage of the aluminum in the stiffener with graphite epoxy always 


results in an increase in o_/P without regard to the percentage 


of composite, size of the oe or the buckling pattern of the 
plate. Thus, it is evident that the torsional stiffness of a rein- 
forced metal stiffener can play a significant role in the stability 
problem and should not be overlooked when high-extensional modulus 


composites are introduced to improve structural efficiency. 


Another popular composite material in developmental use is boron epoxy. 
This composite, although having a higher value of Young's modulus, is 
appreciably more dense and significantly more expensive than graphite 
epoxy. Figure 4 is a plot of the buckling efficiency of'a boron-epoxy- 
augmented stiffened plate. Comparison of Figure 4} with Figure 2 reveals 
that although the use of boron epoxy results in higher values of 2 
than graphite epoxy in the 1,1 buckling configuration, it results in 
greater loss in efficiency when in the 2,2 mode. Higher values of 

SF ./ Diag in the 1,1 mode are the result of the increased Young's 
modulus for boron epoxy, whereas the reduced efficiency in the 2,2 
configuration is attributed to the increased weight of boron epoxy and 


the fact that the shear moduli of the two materials are essentially the 


same. 
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CONCLUDING REMARKS 


This investigation develops a buckling criterion for axially compressed, 
simply supported, flat, rectangular, homogeneous, isotropic plates with 
composite-reinforced, centrally located, longitudinal stiffeners. This 
criterion is obtained from a more general development based on ortho- 
tropic plate theory. Comparison of the results of the present study 
with past investigations shows the derived buckling coefficient to be 
generally more realistic than those which neglect the torsional stiff- 
ness of the stiffener when the plate cross-section aspect ratio between 


stiffeners is at least sixty or less. 


It is also shown that the addition of composites to improve the effi- 
ciency of a stiffened flat plate requires more judicious consideration 
than simply replacing conventional material with high-performance fi- 
brous composites. Effects which may have proved to be negligible in 

the analysis of conventional structures can become significant in the 
conventional-composite combination due to the relatively large spread 
between E and G in advanced composites as compared with that in 
conventional metals. The replacement of a conventional stiffener with 

one consisting of composites will usually lead to lighter structural 
weight, but the undesirable side effects of lower buckling efficiency and 
increased cost may also be introduced. Careful study must be made as to 
the amount of composite to be used, the percentage of increase in buckling 
efficiency actually achieved, and the additional costs involved. It has 
been shown that in certain cases, a rather minimal use of composites will 
result in maximum gain in weight reduction and increase in critical stress 


with relatively small increases in cost. 


When working with advanced composites, the designer of high-performance 
vehicles can no longer neglect the torsional rigidity of a stiffener as 
was normally done with low subsonic-range all-metal systems. Such ne- 
glect can introduce errors which may tend to show that any use of com- 
posites will achieve additional gains in structural efficiency, when in 


fact this is not so. The lower shear modulus of composites requires that 
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the effects of torsional rigidity should always be taken into account, as 
unconservative errors in critical stress may result in either premature 


buckling or reduced buckling efficiency. 


The designer must also be cognizant of all the advantages and disadvan- 
tages of the various composite materials now becoming available. As 
shown herein for graphite and boron epoxy, each of these materials has 
certain advantages, and the most suitable material must be chosen for 
the specific application at hand. Composite materials present the pos- 
sibility of increased strength and stiffness along with decreases in 
structural weight, but careful study must be made with regard to proper 
use and material selectivity to achieve optimized structural efficiency 


and fabricability at minimum cost. 
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APPENDIX I 
DEVELOPMENT OF BUCKLING CRITERION 


The total potential energy, Equation (12) (with Ny = FL ="), upon 
substitution of the deflection function, Equation (15), yields 
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Integration of this expression over the area of the plate yields 
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Now, the first variation of Equation (23) with respect to Aq must 


reduce to zero since any first order change in the total potential 


must vanish for every admissible variation in the lateral displacement 
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Equation (24) reduces to Equation (17) when a factor of mrp. /4L> 





is removed. For the special case of an isotropic square plate 
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Substitution of Equation (25) into Equation (17) yields 
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A second approximation solution of Equation (26) for the one half-wave 


pattern in each direction yields two equations in two unknowns: namely 
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Now, for the Aa S to possess other than trivial values, the deter- 


minant formed by the coefficients of the Bo must vanish. Hence, 
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The solution of this determinant yields Equation (18). The second- 
approximation solution of Equation (26) for the two half-wave pattern 


in each direction yields 
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The solution of Equation (30) yields the result given in Equation (19). 
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APPENDIX II 
DETERMINATION OF EQUIVALENT VALUES OF ele AND GJ 
IN COMPOS ITE-REINFORCED STIFFENER 


9 


Timoshenko and Young~ describe a common method of transforming a com- 
posite beam into an equivalent beam of one material for the determi- 
nation of bending stiffness. This is done by multiplying the width 

of the composite section of the stiffener by the ratio of the Young's 


moduli; that is 


E 
wv 6] See (31) 
1 eo 

base 


where Wy is the original width of the stiffener, and Wo is the 


modified width of the composite section of the equivalent stiffener. 


As an example, suppose EeenD a Der acre 


the equivalent stiffener would consist of one material, the basic 


Then, as shown in Figure 6, 


metal, but it would possess a different shape. The moment of inertia 
of the equivalent stiffener is then found and multiplied by Eades ; 
to give the equivalent value of (EI) for the stiffener. 


, rege nae 
Timoshenko and Goodier : give a formula for the twisting moment or 
torque of a rectangular beam in torsion as 
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However, from Reference 9, the torsional stiffness of a beam in free 


torsion is known to be 


i = (33 ) 


a) 





Substitution of Equation (32) into Equation (443) yields 


Le 6) 
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This value is calculated for each section of the stiffener and multi- 
plied by its respective shear modulus. Addition of these components 


gives the equivalent (GJ) of the stiffener. 
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Figure 1, Geometry of Homogeneous, Isotropic Plate 
With Longitudinal Stiffener. 
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Figure 2. Buckling Efficiency Curve for Aluminum Plate With 
Longitudinal Stiffener Augmented With Varying Per- 
centage (by Volume) of Graphite Epoxy. 
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Figure 43. Buckling Efficiency Curve for Aluminum Plate With 
Longitudinal Stiffener Augmented With Varying Per- 
centage (by Volume) of Graphite Epoxy Ce = Create 
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Figure 4. Buckling Efficiency Curve for Aluminum Plate With 
Longitudinal Stiffener Augmented With Varying Per- 
centage (by Volume) of Boron Epoxy. 
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Geometry of the Transformation of a Composite Beam Into 


an Equivalent Beam of One Material. 
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